Nonequilibrium Steady State in a Quantum System: 
One-dimensional Transverse Ising Model with Energy Current 
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We study the nonequilibrium steady states of an Ising chain in a transverse field, h, by investigating 
the effect of a field, A, which drives the current of energy. The zero-temperature, h — X phase 
diagram is determined exactly and it is found that the energy current appears continuously above a 
threshold, A > Xc{h). The long-range magnetic order (present for h < he, X < Ac) is destroyed in the 
current-carrying state where the correlations are characterized by power-law, oscillatory decay. The 
mechanism which generates the power-law correlations in the current-carrying state is discussed. 
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Nonequilibrium steady states have been much inves- 
tigated but a general theory comparable to the theory 
of equilibrium systems is still lacking. A traditional ap- 
proach to nonequilibrum phenomena starts with the con- 
struction of deterministic phenomenological equations, 
followed by a linear or nonlinear stability analysis of 
the solutions and ends up with the study of the effect 
of fluctuations by adding noise terms to the equations 
10,^ . An alternative and somewhat more microscopic ap- 
proach is the introduction of kinetic Ising models where 
the nonequilibrium steady states are produced by the 
spins being driven by external fields or by being in con- 
tact with several heat baths at different temperatures 



Both of these approaches have lead to successful de- 
scriptions of a number of particular phenomena but they 
have not been very helpful in deducing general conclu- 
sions. The main problem is that the observed proper- 
ties of nonequilibrium steady states are strongly depen- 
dent on the details of the dynamics (even the presence 
or absence of a phase transition or the type of the phase 
transitions, etc. may depend on the details of the tran- 
sition probabilities). From a theoretical point of view, 
the problem is further complicated by the absence of the 
detailed balance condition which provides at least some 
constraints on the possible forms of dynamics of near- 
equilibrium fluctuations. 

A way to avoid the above problem of too much free- 
dom is to investigate systems which, in contrast to models 
such as the Ising model, have natural (quantum mechan- 
ical) dynamics. A nontrivial model with intrinsic dy- 
namics investigated in this paper is the one-dimensional 
(d = 1) transverse Ising model. Once a quantum mechan- 
ical model is chosen, the next question is how to make it 
a nonequilibrium system. A possible route is to couple 



it to heat baths just as in the case of the various kinetic 
Ising models. Unfortunately, the coupling of quantum 
mechanical systems to classical heat baths is not quite 
understood and leads again to much arbitrariness. In or- 
der to avoid this, we follow a different path by noting 
that nonequilibrium steady states are always associated 
with some kind of current (of heat, particle, momentum, 
etc.). Thus we make a nonequilibrium steady state in a 
quantum system by imposing a current on the system. 
For example, in case of the transverse Ising model where 
the only conserved quantity is the energy, we shall con- 
strain the system to have an energy current. It is hard 
to treat, however, a 'microcanonical' type constraint of 
fixed energy fiow and so, we shall actually use the eas- 
ier, 'canonical' description which consists of introducing 
a field which drives the given (energy) current. 

The next question is what to investigate. An obvious 
choice is to study the changes produced by the presence 
of a current in the orderings and in the phase diagram 
of the system. In particular, for the rf = 1 transverse 
Ising model, one can investigate whether the T = Ising 
transition remained in the same universality class and 
whether new phases emerged as a consequence of the en- 
ergy current. Furthermore, it is known that currents in 
nonequilibrium steady states usually generate long-range 
correlations j^]. Thus it is natural to inquire if similar 
phenomena would occur in a quantum nonequilibrium 
steady state. 

The practical realization of the above program consists 
of the following steps: 

i) The starting point is the Hamiltonian which, for the 
d = 1 transverse Ising model, has the form: 
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where the spins are represented by Pauh spin matrices 
cr" (a = X, y, z) at sites £=1,2, of a d = 1 periodic 
chain — erf), and h is the transverse field in units 

of the Ising coupHng. 

ii) The local energy current, Jf, is calculated by taking 
a time derivative of the energy density, using the quan- 
tum mechanical equation of motion, and representing the 
result as a divergence of the energy current (fi. = 1 is 
used): 

Ji = \<T>u-<^hi) ■ (2) 

iii) The 'macroscopic' current J — Ji is added to 
Hi with a Lagrange multiplier, —A, 

H = Hi-XJ . (3) 

Note that the energy current, J, is associated with Hi 
and not with the new Hamiltonian, H . We also empha- 
size that H is just another equilibrium Hamiltonian, it 
differs from Hi by an extra term which breaks the left- 
right symmetry of Hi. Finding the ground state of H, 
however, gives us the minimum energy state of Hi which 
carries an energy current, J = {J/N) (brackets () de- 
note the expectation value in the ground-state of H). 
Thus the ground-state properties of H provide us with 
the properties of the nonequilibrium steady states of the 
transverse Ising model. 

iv) The correlations and ordering properties of the sys- 
tem are investigated as a function of A. Of course, one 
can also find that value of A which generates a given en- 
ergy current, J = (J/N), in the transverse Ising model 
and then all steady state properties can be obtained in 
terms of J. 

It should be noted that the terms in J can be re- 
arranged so that the energy current becomes a current 
through the bonds: 

^ = E-^^ = iE(^?^'+i-'^K+i) • (4) 

e e 

This form is remarkable because Ji can be recognized 
as the Dzyaloshinskii-Moriya interaction in the theory of 
weak ferromagnetism Moreover, the current of the 
z component of the magnetization in the transverse XY 
model is also of the above form 0, and similar terms 
(apart from an overall factor of i which makes the Hamil- 
tonian nonhcrmitian) appear in the description of driven 
diffusive systems as well as in the Hamiltonian de- 
scription of the direct electric field of the 6- vertex model 

We turn now to the solution of the problem. Since 
[Hi, J] — 0, the Hamiltonian H is diagonalized using the 
same transformations which diagonalize if/ [|lO|,^: first, 
creation-annihilation operators are introduced (a^ = 
cr^ ± ia^), then the Jordan- Wigner transformation ||ic| ] 
is used to transform them into fermion operators (q, c^) 



and, finally, a Bogoljubov transformation jlj] is em- 
ployed on the q and —q components of the Fourier trans- 
forms of Q-s. As a result one finds the spectrum of exci- 
tation energies as tjg = |Aq| where 

^9 = ^ (Vl + + 2hcosq + (ainq^ . (5) 

Here the wave numbers are restricted to — tt < g < tt in 
the thermodynamic limit {N — > oo) and we introduced 
C = Xh/2 which appears to be the natural variable in- 
stead of A. Fig.|l| displays the spectrum for h = 0.5 and 
various ^ and one can see that the q —q symmetry 
of the spectrum is broken for C 0- For small C, when 
Ag > 0, however, the ground-state remains that of the 
transverse Ising model (C = 0) and, accordingly, no en- 
ergy current flows (J = 0). This rigidity of the ground 
state against the symmetry-breaking field which drives 
the energy current is a consequence of the facts that the 
fermionic spectrum of the transverse Ising model has a 
gap and that the operator J commutes with Hi. Actu- 
ally, the rigidity persists even at the critical point he = 1 
where the gap disappears. The reason for this is more 
subtle and related to the fact that the critical spectrum 
of Hi has a finite slope at g = ±7r. 
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FIG. 1. Spectrum of the transverse Ising model in the pres- 
ence of a field (A) which drives the current of energy. The 
excitation energies are given as oj, = |Ag|. The qualitative 
picture is the same at all transverse fields h except at ft = 1 
where the q — ±7r gap disappears. 



The ground-state properties do change when A^ < in 
an interval [q- , 17+] since these q states become occupied 
and q± depends on both h and C. The line Cc(^) which 
borders the region of unchanged transverse Ising behavior 
is obtained from the conditions 

A, = , dA,/dq^O , (6) 

and the solution 
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h > 1 
/i < 1 



(7) 



is displayed on the phase diagram (FigJ|) as a solid line. 
Fig.^ also shows the phase boundary (dashed line) be- 
tween the magnetically ordered {h < 1, ( < 1) and dis- 
ordered (/i > 1, < /i) transverse Ising regions. Since 
the ground state is independent of C for C < Cc, the tran- 
sition across the dashed line is a second order transition 
belonging to the d = 2 Ising universality class 

In the following we shall argue that the C > Cc region 
can be considered as a distinct phase since i) the energy 
current is nonzero, ii) there is no long-range magnetic 
order, and iii) the magnetic correlations are oscillatory 
with amplitudes decaying as a power of distance. 
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FIG. 2. Phase diagram in the h — ( plane where h is 
the transverse field while ^ — hX/2 is the effective field 
which drives the current of energy. Power-law correlations 
are present in the current-carrying phase (J ^ 0) and on the 
Ising critical line (dashed line). 



The energy current is expressed as a product of two 
fermion operators and can be calculated exactly. It is 
zero for C ^ Cc while, for C > Cc, one obtains: 



J = (J/N) 
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At the phase boundary, Cc we have J ^ as 

J ' 



C-Cc 



h 1 

h = 1 



(8) 



(9) 



Thus the Cc line can be considered as a line of second 
order transition provided J is viewed as an order param- 
eter. 

There are several other quantities which change non- 
analytically as we cross into the current carrying phase. 



For example, the change in the z component of the mag- 
netization at fixed h is given by 



(C-Cc)^/^ 

C-Cc 

(C-Cc)^/^ 
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h = 1 
h < 1 



(10) 



The different exponents obtained for h > I and h < 1 
reflect the fact that the J ^ phase neighbors magnet- 
ically ordered ((erf) ^ 0) and magnetically disordered 
((erf) = 0) phases along the C^^" and C^T portions of the Cc 
phase boundary. 

The long-range magnetic order ((erf) ^ 0) disappears 
for C > Cc- This can be seen by investigating the {crfcrf+n) 
correlations which can be expressed through Pfaffians 
[|l3|,^. In the presence of long-range order, one should 
have (erf CTf^„) — > (erf )^ ^ for rt ^ oo. Instead, we find 
that the correlations decay to zero at large distances as 



\/n 



■ cos{kn) 



(11) 



where the wavenumber, k, depends only on C, 

k ~ arccosC^^ . (12) 



The above results (11 12) are exact in the ( oo limit 
where Q{h, oo) = e^T^T^/^A'^ « 0.147 is the amplitude 
of the {cr^cr^^,^) correlations of the d — 1 XY model |l^] 
{A « 1.282 is Glaisher's constant). The connection to 
the XY model and the details of the calculation will be 
published separately |15||. 
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FIG. 3. The correlation function (erf (jf^„) at various trans- 
verse fields h with C, = h\/2 fixed at ^ = 2. The lines are fits 
to the large n behavior which is assumed to be of the form 
Qn~^^^ cos kn. Changes in the amplitude Q occur mainly 
near the phase boundary {h — > 2) where Q ^ 0. Note that 
the convergence to the asymptotic behavior becomes slower 
as ft — + 2. For comparison, we have also plotted the correla- 
tions {o"f (T^^„) ~ n"^''^ on the Ising critical line (dashed line 
on Fig.2). 
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For < C < oo, the expressions ( |Tl| , p^ ) are obtained 
from the numerical analysis of (cJcr^Ln) for finite n < 
100. Typical results are shown on Fig.|^. The numerical 
results are behaved well enough so that, in addition to 
the n dependence of {creo-i+n)^ functional form of 
the amplitude Q{h, () can also be investigated. We find 
that the following expression fits the results everywhere 
(C > Cc) except close to the phase boundaries, (f, or 
near the h — I line: 

/^2_;j2X 1/4 

Q(/i,C)«Q(/i,^) ■ 

Near either ( — or /i = 1, the convergence to the 
asymptotic form (O) is slow and the numerics becomes 
unreliable. The convergence becomes fast again on the 
h — 1 line itself where the following form appears to be 
valid: 

1/8 

g(i,C)«g(^^)(^^5^) . (14) 

Power-law decay of correlations for C > is present 
in other physical quantities as well. For example, the en- 
velopes of both (cr|CT|_,_„) and {JiJi+n) decay in the large 
n limit as |p^ . Thus we find power-law correlations 
in the current carrying state in agreement with the no- 
tion that power-law correlations are a ubiquitous feature 
of nonequilibrium steady states . 

Power-law correlations in quantum models are usually 
associated with a gapless excitation spectrum. Thus, pro- 
vided the transverse Ising model can be considered as an 
instructive example, we can see a general connection be- 
tween the emergence of power-law correlations and the 
presence of a current. Indeed, let us assume that a sys- 
tem with Hamiltonian has a spectrum with a gap 
between the ground-state and the lowest excited state. 
Furthermore, let J be a 'macroscopic' current of a con- 
served quantity such that [Ho, J] — 0. Generally, there is 
no current in the the ground state and adding — AJ to i?Q 
does not change the ( J) = result for small A. Current 
can flow only if some excited states mix with the ground 
state and, consequently, a branch of the excitation spec- 
trum must come down and intersect the ground-state en- 
ergy in order to have (J) ^ 0. Once this happens, how- 
ever, the gap disappears and one can expect power-law 
correlations in the current-carrying state. 

Admittedly, the above argument is not strict and is just 
a reformulation (in general terms) of the transverse Ising 
results. We believe, however, that the above picture is 
robust and suggestive enough to warrant further studies 
in this direction. 

In summary, we have seen that nonequilibrium steady 
states in quantum systems can be constructed by switch- 
ing on fields which drive the current of conserved quanti- 
ties. Studies of these states may be of importance since 
their construction is free from the arbitraryness inherent 



in classical nonequilibrium dynamics. On the particu- 
lar example of the transverse Ising model, we could ob- 
serve that there are power-law correlations in the current- 
carrying states and we argued that this may be a general 
feature of quantum systems. 
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